Butterfly hysteresis loop at non-zero bias field in antiferromagnetic molecular rings: 

cooling by adiabatic magnetization 
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At low temperatures, the magnetization of the molecular ferric wheel NaFe6 exhibits a step at a 
critical field B c due to a field-induced level-crossing. By means of high-field torque magnetometry 
we observed a hysteretic behavior at the level-crossing with a characteristic butterfly shape which is 
analyzed in terms of a dissipative two-level model. Several unusual features were found. The non- 
zero bias field of the level-crossing suggests the possibility of cooling by adiabatic magnetization. 



One of the most exciting recent discoveries in mag- 
netism is that of magnetic hysteresis of purely molecular 
origin m metal complexes like Mn i2 or Fe 8 @ |, |, !,§,§■ 
In these large-spin molecules spin reversal is blocked at 
low temperatures due to an energy barrier. Recently, a 
different type of molecular hysteresis has been observed 
in the molecule called V15 0. Here, the spin system 
is dynamically driven out of equilibrium leading to hys- 
teresis loops with a characteristic butterfly shape. The 
relevant time scale is on the order of seconds because of 
the phonon bottleneck effect. All in all, V15 is an experi- 
mental realization of a two-level system with dissipation 

0- 

In this work, we investigated the molecular ferric wheel 
[NaFe 6 L 6 lCl-6CHCl 3 with L = N(CH 2 CH 2 0)3, or NaFe 6 
in short [8|. Ferric wheels attracted much interest be- 
cause of the possibility of coherent quantum tunneling 
H . In NaFe6, the six iron(III) ions form a ring of spin- 
5/2 centers with dominant antiferromagnetic nearest- 
neighbor Heisenberg interactions (ll). At zero field the 
ground state is nonmagnetic (S — 0), but application 
of a magnetic field leads to a level-crossing at a critical 
field B c where the ground states changes abruptly to the 
first excited S = 1, M = — 1 level [inset of Fig. 1(a)] 
p2| , Thus, from a general point of view, NaFe6 rep- 
resents a two- level system for fields close to B c . We will 
show that, for appropriate conditions, NaFeg exhibits a 
butterfly hysteresis similar to that observed in V15. This 
allows for important conclusions concerning fundamen- 
tal physical aspects of ferric wheels. Furthermore, we 
found several striking features of the spin reversal pro- 
cess, demonstrating that the dissipative two-level system 
realized in NaFe6 is actually a very distinctive one. 

The magnetic properties of single crystals of NaFeg 
were investigated by high-field torque-magnetometry. 
The single crystals were prepared as reported f|. They 
crystallize in the space group R3. The cation [NaFe6L@] + 
exhibits Sq molecular symmetry with the Se symmetry 
axis oriented perpendicular to the plane of the cluster. 
The complex shows strict uniaxial magnetic behavior 
pi The uniaxial axis coincides with the molecu- 

lar Sq symmetry axis. The torque was detected with 



homemade cantilever devices micromachined from crys- 
talline silicon Resolution was typically 10~ n Nm. 
For measurement, a single crystal was selected by light 
microscopy in the mother liquor, was directly put from 
the solution into Apiezon grease and mounted on the 
cantilever. The weight of the crystals investigated was 
typically 20 fig, non-linearity was less than 1%. The 
torquemeter was inserted into a 17 T cryomagnet with a 




FIG. 1: Field dependence of the torque of a NaFe6 single crys- 
tal demonstrating the butterfly hysteresis loop as it typically 
appears for field sweep rates in the saturation regime (see 
text). The inset in panel (a) schematically shows the field de- 
pendence of the low-lying energy states and the level-crossing 
at B = B c leading to the pronounced torque step. In panel 
(b), the dashed curve represents the measured equilibrium 
curve, the inset sketches the behavior of the spin tempera- 
ture T s when the field is swept horn B > B c . 
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FIG. 2: Measured hysteresis loops for (a) three field sweep 
rates at T = 0.32 K and (b) three temperatures at a constant 
sweep rate of 23 mT/s. Only the branch with B > B c is 
plotted. 

vacuum loading 3 He insert. The angle between the crys- 
tal's uniaxial axis and magnetic field could be aligned in 
situ with an accuracy of 0.3°. Calibration of the torque 
signal is accurate to ±15%. Temperature was measured 
with a RuC>2 thick film resistor to within ±30 mK. 

A typical result for the field dependence of the torque 
at T — 0.32 K is shown in Fig. 1. The negative sign of the 
torque manifests the hard axis anisotropy of NaFe@ [|ll| . 
The torque signal exhibits a steplike field dependence at 
B c due to the abrupt switch of the ground state from the 
zero-field S = ground state to the first excited 5* = 1 
state [inset of Fig. 1 (a)] . This phenomenon is well known 
for ferric wheels jyj |lj, |l3| , but for the first time we also 
observed a pronounced hysteresis of the torque signal for 
fields close to B c , as shown in more detail in Fig. 1(b). 
The dependence of the butterfly curves on the sweep rate 
of the field and on temperature is presented in Fig. 2. The 
field position of the curves, i.e. the value of B c , showed 
the expected angular dependence [p"T|| , but otherwise we 
didn't observe an effect of the angle. 

The underlying physics is quickly identified. Mea- 
surements performed with torque sensor and crystal im- 
mersed in 3 He exchange gas resulted in perfectly re- 
versible torque curves. In this work, however, we used a 
vacuum loading 3 He insert and thermal contact between 
crystal and bath was established via the silicon cantilever. 



In this configuration, thermal coupling should be weaker 
than with exchange gas. Obviously, the hysteretic behav- 
ior shown in Fig. 1 reflects a non-equilibrium state due to 
limitations in heat flow indicating that phonons play an 
important role in the process of changing the spin state 
at the level-crossing. 

Similar butterfly hysteresis curves were observed re- 
cently for the molecular cluster V15 and were success- 
fully described in terms of a phonon bottleneck scenario 
(PBS) @, Q. The essence of the PBS is that energy ex- 
change between spin system and bath occurs via phonons 
fL5|| , The phonon bottleneck arises because the number 
of spins is much larger than the number of available res- 
onant phonon modes at low temperatures. That is, the 
heat capacity of the spins C s is much larger than that of 
the phonons C p (b = C s /C p ^> 1). Therefore, the phonon 
temperature T p follows that of the spin system T s very 
quickly (on times t\ xs b~ 1 t s ) while energy transfer from 
the phonons to the bath is drastically delayed (time scale 
t 2 ~ btp). Here, t s and t p denote the spin-phonon and 
phonon-bath relaxation times, respectively. For the time 
scales of the present experiments, the spins and phonons 
may be regarded as a single coupled system (T s = T p ) 
which is weakly coupled to the bath (T s ^ T). 

The butterfly hysteresis curves reflect the behavior of 
the spin temperature in the PBS. Spin temperature and 
torque are monotonously related: for a given magnetic 
field, I -7- 1 is the larger the smaller T s if B > B Cl and 
vice versa if B < B c . Approaching the level-crossing 
from a starting field of e.g. B > B c , the spin temper- 
ature first decreases since the spin-phonon system can- 
not pick up heat quickly enough from the bath [inset of 
Fig. 1(b)]. Thus, for sufficiently fast sweeping rates, the 
spin-phonon system is virtually decoupled from the bath 
resulting in a (quasi) adiabatic demagnetization process. 
At B w B Cl the spin system reaches a minimal tem- 
perature T ' s ,min- With further decreasing field, the spin 
system heats up (reverse of adiabatic demagnetization) 
but overshoots the bath temperature since it now receives 
heat from the bath which started to flow to reestablish 
equilibrium but reaches the spin system delayed due to 
the finite relaxation time ti. 

The phonon bottleneck effect is obviously the under- 
lying physical mechanism for the observed hysteresis in 
NaFe6, and in so far the situation is analogous to that of a 
conventional two-level system. However, there are some 
particular characteristics to be noted, some of which are 
quite striking: 

(1) With our setup we observed "heat-flow" hysteresis 
effects only for NaFe6 (and related ferric wheels). We 
have investigated quite a number of different molecular 
spin clusters, but for them the thermal coupling provided 
by our silicon sensors was sufficient to achieve reversible 
torque curves. In this respect, the observation of hystere- 
sis in NaFe6 is notable. 

(2) A distinctive feature of the two-level system real- 
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FIG. 3: Dependence of the minimal spin temperature T s , m in 
on the field sweep rate for two samples. Dashed lines are 
guides to the eye. The inset schematically presents the two 
states at the level-crossing and indicates the splitting Ao dis- 
cussed in the text. 



ized in NaFeg is that the two states involved belong to 
different spin and spatial quantum numbers. In NaFe@ 
the strong exchange limit is realized and the total spin 
quantum number S is a good (though not exact) quan- 
tum number [juj Additionally, the cyclic molec- 
ular symmetry implies a spin permutational symmetry 
[16 1 and states may be classified by the wave vector k 
[13, |l7| ]. The two states relevant in this work belong to 
S = 0, k = it and S = 1, k = 0, respectively, so that 
phonon absorption/emission is constrained by A5 = 1 
and | Afc| = tt. To the best of our knowledge, the effects of 
these restrictions for the energy exchange between spins 
and bath have not been investigated so far. But con- 
ceivably the number of available phonons may be further 
reduced by these constraints. The Ak selection rule may 
be of particular relevance since it holds exactly as long 
as the molecular symmetry exhibits a C*2 symmetry axis. 

(3) With a numerical analysis of the hysteresis curves 
based on the PBS, Chiorescu et al. evidenced a zero- field 
splitting for the molecule V15 0. In the case of ferric 
wheels, within a description in terms of a spin Hamilto- 
nian consisting of terms related to the electronic spins 
only, the level-crossing at B = B c is a true crossing by 
virtue of translational symmetry ^|. However, inter- 
actions with e.g. the nuclear spins may lead to a splitting 
of the levels at the crossing point |l^, |l9| . Their identifi- 
cation would be of considerable importance for the issue 
of quantum coherence in antifcrromagnetic rings since 
they may act as sources of decoherence [jOl. But so far, 
no clear evidence for a splitting of the level-crossings, i.e. 
for anti-crossings, has been reported. 

We thus analyzed our data with the numerical proce- 
dure developed in Ref. as adapted to torque mea- 
surements. Although we could generate curves which 
roughly resembled the experimental ones, we could not 
reproduce the data for different temperatures and sweep- 
ing rates satisfactorily. Values for e.g. Aq, the splitting 



of the level-crossing at B = B c (inset of Fig. 3), showed 
unphysical scattering. This indicates, that the PBS as 
outlined above is a too simple picture. 

Therefore, we considered the dependence of the min- 
imal spin temperature T s ^ m i n as function of the sweep 
rate of the magnetic field. In conjunction with a model 
Hamiltonian, T s>m i n could be determined reliably from 
the slope of the torque curves at B = B c as the curves 
are essentially reversible in this field regime. In this way, 
T s ,min could be estimated with a reproducibility of 5 mK, 
but the values are subject to a systematic error of about 
±15 mK due to inaccuracies in the calibration of the 
torque signal. The results for two samples are shown 
in Fig. 3. For sample A, T Sy7nin first decreases linearly 
with increasing sweep rate and then saturates at about 
0.17 K. It should be noted, that the value of T s ^ nin — 

0. 325 K extrapolated to zero sweep rate agrees well with 
the bath temperature as measured with the RuC>2 sen- 
sor. In the saturation regime, the spin system obviously 
behaves quasi adiabatically and it is thus very suggestive 
to identify the saturation temperature as a level-splitting, 

1. e. A = 0.17 K. 

We have also analyzed the broadening of the torque 
step as measured under reversible conditions for various 
temperatures down to 0.4 K. This method provides an 
alternative way to estimate Aq [O, 18 1. However, even 



at the lowest temperatures the width of the step is in full 
accord with thermal broadening resulting in an upper 
limit of 0.05 K for A . 

We have no explanation for this serious discrepancy. It 
suggests that the minimal spin temperature in the limit 
of infinite sweep rates exceeds the thermodynamically de- 
termined splitting of the level-crossing. The PBS in the 
present form obviously does not account for all the ef- 
fects responsible for limiting T Stm i n and/or broadening 
the torque step. It should be noted that the discrepancy 
cannot be explained by "extrinsic" effects as e.g. crystal 
imperfections (twinning, mosaicity, etc.) |is}| . 

Incidentally, Fig. 3 demonstrates the role of the ther- 
mal coupling between crystal and bath. In case of sample 
B it is obviously better than for sample A as the crossover 
is shifted to larger sweep rates. 

(4) As mentioned already, T a>m i n as function of the 
sweep rate exhibits a crossover from a linear regime to 
a saturation regime (Fig. 3). This is readily under- 
stood within the PBS. However, in the case of NaFeg 
the crossover is accompanied by a striking change of the 
hysteretic behavior. In the saturation regime the curves 
are symmetrical around B c [Fig. 1(b)], but in the linear 
regime the hysteresis becomes clearly asymmetric. Two 
example curves are shown in Fig. 4. The asymmetry is 
more pronounced the slower the sweep rate or the higher 
the (bath) temperature and was always found to be such 
that the hysteresis is larger for B > B c than for B < B c . 

It seems as if energy exchange between spin system and 
bath is different for B < B c and B > B c . This is not eas- 
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FIG. 4: Two example curves showing the asymmetry of the 
hysteresis loop around B c discussed in text. 



wheels. The present work suggests a novel experimental 
route to measure them. 

For slow sweep rates the hysteresis loop exhibited a 
striking asymmetry which is unexplainable within a stan- 
dard phonon bottleneck scenario. One may only specu- 
late about the role played by the unusual quantum num- 
bers of the states involved in the level-crossing. It also 
has been demonstrated that the NaFeg system allows for 
cooling by adiabatic magnetization, i.e. for the opposite 
of what is generally known. At the level-crossing the fer- 
ric wheel NaFe6 represents a dissipative two-level system 
- but with quite unusual properties sheding new light on 
a model of which one thought to know all. 



ily explained. It is hardly conceivable that the thermal 
coupling between phonons and bath depends on magnetic 
field. But then one has to conclude that the asymmetry 
is related to peculiarities of the spin-phonon coupling in 
NaFe@. However, the arguments presented above and in 
Refs. @ aim to show that spin-phonon coupling is 
so fast that in our experiments T s = T p is always main- 
tained so that any details of the spin-phonon coupling 
are indiscernible! As demonstrated by the symmetry of 
the curves for large sweep rates and under reversible con- 
ditions [Fig. 1(b)], an asymmetry of the energy levels is 
not an issue. 

(5) The most obvious distinctive feature of the NaFeg 
two-level system is the non-zero bias field of the cross- 
ing, i.e. of B c > 0. As shown in the above, the spin 
system cools down by approaching the critical field B c as 
it is quasi adiabatically decoupled from the environment 
for large enough sweep rates. If one starts at B > B c 
and decreases field, the situation corresponds to cooling 
by adiabatic demagnetization. However, if one starts at 
B < B c (e.g. B = 0) and increases field, the reverse 
situation is realized, namely that of cooling by adiabatic 
magnetization. Of course, in the present work no real adi- 
abatic conditions were realized - but the effect is clearly 
evident from Fig. 1. 

In summary, we observed a butterfly hysteresis loop 
of the torque at the first ground state level-crossing in 
the molecular ferric wheel NaFe6- The results show that 
phonons are mostly responsible for transitions between 
the two states near the level-crossing. This clarifies a re- 
cent proposal where it has been argued (as hysteresis was 
not yet observed) that dipole-dipole coupling between nu- 
clear and electronic spins plays an important role [p0[ . 
The minimal spin temperature T Sj7n i„ achieved during 
the hysteresis loop was found to saturate for large sweep 
rates of the magnetic field - the first direct experimen- 
tal evidence for interactions of the electronic spins with 
the environment in ferric wheels. A detailed analysis of 
the processes determining T Stm i n should provide deep in- 
sight into the conditions for quantum coherence in ferric 
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